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Abstract 

In recent lattice QCD studies, the Feynman-Hellmann theorem is often used to estimate separate 
contributions of the connected and disconnected diagrams to the nucleon sigma term. We demon- 
strate through a simple analysis within an effective model of QCD why this could be dangerous 
although the theorem is naturally expected to hold for the sum of the two contributions, i.e. the 
net nucleon sigma term. 
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I. INTRODUCTION 



The nucleoli sigma term T^t^n is believed to be a quantity of fundamental importance 
in that it gives a measure of the explicit chiral symmetry breaking of QCD. In fact, it 
characterizes the effect of finite quark mass on the mass of the nucleon as M^v = Mq + E^^r , 
with Mo being the nucleon mass in the chiral limit. Recently, the JLQCD collaboration 
reported an estimate of the nucleon sigma term with good precision based on the the overlap 
'ermion action, which preserves exact chiral symmetry and flavor symmetries on the lattice 
l|. They estimated the separate contributions of the connected and disconnected diagrams 
to the nucleon sigma term by utilizing the Feynman-Hellmann theorem derived within the 
framework of the partially quenched QCD (PQQCD), where the quarks that couple to 
external sources for the asymptotic hadrons, i.e. the valence quarks, are distinguished 
from those that contribute to the quark determinant, i.e. the sea quarks. They found 
that the connected diagram gives a dominant contribution to the nucleon sigma term and 
the disconnected-diagram contribution to it is fairly small. It appears to contradict our 
experience within the chiral quark soliton model (CQSM), in which we found the dominance 
of the Dirac-sea quarks over the valence quarks in this special observable {2] (From 
the physical ground, the valence and Dirac-sea contributions in the CQSM is expected to 
correspond to the connected- and disconnected-diagram contributions in the lattice QCD, 
at least approximately.) 

What is the cause of this discrepancy ? There appears to be httle reason to suspect the 
validity of the Feynman-Hellmann theorem, especially because it can be proved on quite 
general theoretical postulates. At the same time, however, one should recognize the fact 
that the general proof of the theorem in textbooks of quantum mechanics is given only 
for the total mass or the total Hamiltonian. (See jsj], for instance.) If one divides the 
total contribution into two parts, it is highly nontrivial whether the theorem holds for the 
individual pieces separately. One might think that it is not a serious problem, since the 
Feynman-Hellmann theorem is anyhow expected to hold for the net nucleon sigma term and 
since only the sum is a quantity of physical interest. However, the authors of [l| made a semi- 
quenched estimate of the strange quark content of the nucleon within the same framework of 
two-flavor QCD utilizing the Feynman-Hellmann theorem, thereby being led to a remarkable 
conclusion that the ss components in the nucleon is very small in contrast with several past 
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estimates in the lattice QCD [1Q|]-[12|. Whether this estimate is justified or not may depend 
on whether the use of the Feynman-Hellmann theorem for separating the connected- and 
disconnected-diagram contributions to the nucleon sigma term is justified or not. 

The purpose of the present paper is to show why a naive apphcation of the Feynman- 
Hellmann theorem can be dangerous when it is used for the separation of the nucleon sigma 
term into the two pieces. The strategy for verifying our claim is as follows. First, we recall 
the fact that, within the framework of the CQSM, we can directly calculate the separate con- 
tributions of the valence and Dirac-sea quarks to the nucleon sigma term, thereby confirming 
that the latter is dominant over the former. Second, we show that a naive application of the 
Feynman-Hellmann theorem leads to a totally different answer from the direct calculation, 
although the sum of the valence and Dirac-sea contributions are exactly the same in the 
two ways of calculating the nucleon sigma term. Next, we shall show that careful inspection 
of the derivation of the theorem indicates the necessity of a correction term, which fills up 
the gap between the direct calculation and the naive application of the Feynman-Hellmann 
theorem. Finally, bearing in mind our finding in the CQSM analysis, we shall make some 
remarks on the corresponding analysis of the nucleon sigma term in the lattice QCD by 
using the Feynman-Hellmann theorem. 

II. THE NUCLEON SIGMA TERM IN THE CHIRAL QUARK SOLITON MODEL 
A. A direct calculation 

We begin with the effective Lagrangian of the chiral quark soliton model (CQSM) with 
an explicit chiral symmetry breaking [2|,[3| : 



(1) 



where Cq denotes the chiral symmetric part given by 




(2) 



with M being the dynamically generated quark mass, and 




(3) 
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On the other hand, 

£' = — moi/j{x) ipi^x), (4) 

is thought to simulate a small deviation from the chiral symmetry limit with mo being the 
bare quark mass. Note that the effective quark mass in the physical vacuum ([/ = 1) is given 
by M = M + rriQ. The model contains the parameters, M, mo, fn and some physical cutoffs. 
Throughout the present study, we set f-j^ = 93 MeV and M = 375 MeV, while the bare 
quark mass mo is varied around the reference value mo = 6 MeV. To fix the regularization 
parameters, we first define the effective action Seff[U] through the relation 

Z = J VnV^V^P^ e'l '^'''^ = J Ptt e'^^^t^l (5) 

Next, to get rid of ultraviolet divergences contained in this definition, we introduce the 
regularized effective action in the proper-time regularization scheme by 

SITAU] = I [ ^ ^(r) Sp (e'^^^^ - e--^^^«) , (6) 

with 

D = i ^ - MW^ - mo, Do = id - (M + mo). (7) 

The regularization function (f{T) is introduced so as to cut off divergences appearing as a 
singularity at r = 0. For determining it, we require that the regularized theory reproduce 
the correct normalization of the pion kinetic term as well as the mass term. Using the 
standard derivative-expansion technique, this gives two conditions : 



4:7T^ Jo T 

NcM f-- dr ^^^2 2 (r,^ 

-,v{r)e =m^. (9) 

Since Schwinger's original choice v^(t) = 6{t — 1/A^), with A being a physical cutoff energy, 
cannot satisfy the above two conditions simultaneously, we use a slightly more complicated 
form as 

(^(r) = ce[r-l/K\) + (l-c)0(r-l/A2) (10) 



with c = 0.720, Ai = 412.79 MeV and A2 = 1330.60 MeV. The soliton construction in the 
CQSM starts with a static mean- field configuration of hedgehog shape as fl^.fisl 

U^o'{x) = e'^'''^-'^^^'\ (11) 



The quark field in this mean-field obeys the Dirac equation : 



with 



H\n) = En\n), 



H = ^ + M/3e^^5'^-^^('-) + mo/5. 



(12) 



(13) 
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FIG. 1: The schematic energy spectra of the Dirac equation with the hedgehog mean field {x) 



A characteristic feature of this Dirac equation is that one deep (single-quark) bound state 
appears from the positive-energy Dirac continuum. We call it the valence quark orbital. An 
object with baryon number one with respect to the physical vacuum is obtained by putting 
Nc{— 3) quarks into this valence orbital as well as all the negative-energy (Dirac-sea) 
orbitals. Accordingly, the total energy of this baryon-number-one system is given by 



Estatic\U\ — Eyal\U] + Esea\U]. 



Here Eyai represents the valence quark contribution to the static energy, i.e. 



EUU] = N,E^[Ul 



(14) 



(15) 



with Eq being the eigen-energy of the valence quark level. On the other hand, Egea stands 
for the energy of the polarized Dirac sea. Regularizing it in the proper-time scheme, we have 



(16) 



The energy of the physical vacuum (f/ = 1) is subtracted here with being the eigen-energy 
of the vacuum Hamiltonian Hq = H[U 1]. The most probable pion-field configuration (or 
the self-consistent mean field) is determined by requiring the stationary condition for the 
total energy, 

j^^E,,^uc[Fir)] = 0. (17) 
This Hartree problem with infinitely many Dirac-sea orbital can be solved by using the 



numerical method of Kahana, Ripka and Soni [16|,[17|. It also enables us to evaluate any 
nucleon observables with full inclusion of the Dirac-sea quarks. Of our particular interest 
here is the nucleon sigma term, which is defined as a J2tvN = "^o ^ with rriQ being the current 
quark mass and a being the nucleon scalar charge given as a = {Nluu + dd\N). By taking 
care of the consistency with the basic equation of motion of the model, the regularized 
expression for the nucleon scalar charge a is given as 

O' = d-^al + ^sea, (18) 

where 

^val = iVc(0|a°|0), (19) 

asea = — :r X! -Fi^ri) {n\a'^\n) — (vacuum subtraction), (20) 
^ n 

with the regularization function, 

1 rOO rjj- „ 

HEn) = ^ ^y.(r)E.e-^"^ (21) 

Numerically, we find that 

a ~ 6.86, (22) 

with 

a,al ^ 1.91, (Tsea ^ 4.95, (23) 

which clearly show the dominance of the Dirac-sea contribution over the valence quark one. 
(With the choice mg = 6MeV, the above nucleon scalar charge gives J^-kN — 41.2 MeV. 
We recall that the nucleon scalar charge, especially its Dirac-sea contribution, is a quantity 
which is extremely sensitive to the regularization scheme. The Pauli-Villars regularization 
scheme, which is also used frequently in the CQSM, leads to much larger nucleon sigma term 
ranging from 48 MeV to 72 MeV depending on the bare quark mass mo 4], 8|.) Anyhow, 



the predictions of the CQSM shown above appears to be quahtatively consistent with the 
results of the older simulations in quenched lattice QCD by based on the Wilson quark action 
[l^l'Hl- (These old calculations by using the Wilson-type fermion, which violate the chiral 
symmetry on the lattice, were criticized, however. The criticism is that such calculations 
can give rise to a significant lattice artifacts in the sea quark content arising from the sea 
quark mass dependence of the additive mass renormalization and lattice spacing [l9|,|l|.) It 
however appears to contradict the recent results of JLQCD collaborations by utilizing the 
Feynman-Hellmann theorem within the framework of the overlap fermion, which indicates 
the dominance of the contribution of the connected diagram over that of the disconnected 
one l|]. What is the cause of this discrepancy? To answer this question, we think it useful 
to evaluate the nucleon scalar charge by utilizing the Feynman-Hellmann theorem within 
the same CQSM. 

B. A naive application of the Feynman-Hellmann theorem 

We begin with the general statement of the Feynman-Hellmann theorem. The theorem 
states that 

Ie ^ (24) 

where 

• H{a) is a Hamiltonian operator depending on a continuous parameter a. 

• I "^{a)) is an eigenstate of the Hamiltonian, depending implicitly upon a. 

• E is the eigen-energy of the Hamiltonian H{a). 

In our present application, the bare quark mass rriQ plays the role of the parameter a, 
and the Hamiltonian is given by 

H{mo) = + M/3e'^»'^'-^^(^) + mo(3. (25) 

Thus, we obtain 

1^ - ^ - ^° P'') 

On the other hand, the eigenstate are given as 

l^("^o)) = n 4 \vac), (27) 

n&occ 



where represents the creation operator that creates a quark in the single-quark eigen- 
state \n) of the Hamiltonian H{mQ), while \vac) is the corresponding empty vacuum. The 
Feynman-Hellmann theorem then dictates that 

^^E{mo) = (^(mo)|7°l^(^o)). (28) 

Since the r.h.s. is nothing but the scalar charge a of the nucleon, we immediately get 

(J = T^^(mo), (29) 

which is the anticipated result. Remembering that the total energy is given as the sum of 
the energy of the valence quarks and that of the Dirac-sea quarks, one would further expect 
that 

O' = a^al + ^sea, (30) 

with 

d 

O-val = Kalimo), (31) 

d 

(Tsea = 7^ -Esea("^o)- (32) 

Within the CQSM, we can solve the eigenvalue problem for any value of mg to obtain E^ai 

and Esea as functions of mo, so that we can readily calculate the r.h.s. of ( l3T|l and (l32ll . 
(The mean field, or the soliton profile function F{r), is fixed throughout this calculation.) 
In that way, we obtain 

a ~ 6.87, (33) 

with 

a^ai ^ 11.18, asea = -4.31. (34) 

This should be compared with the answer of the direct calculation of the nucleon scalar 
charge described in the previous subsection : 

^^""^ = ^iai + ^iel ^ 1-91 + 4.95 ^ 6.86. (35) 

One finds that the two ways of calculating the nucleon scalar charge give totally different 
answers for the individual contributions of the valence and Dirac-sea quarks. Nevertheless, 
both give practically the same answer for the sum of the two contributions, i.e. for the 
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net scalar charge of the nucleon, or equivalently for the net nucleon sigma term. Roughly 
speaking, the valence and Dirac-sea contributions in the CQSM corresponds to the connected 
and disconnected-diagram contributions in the lattice QCD. Then, what should be clarified is 
the reason why the naive application of the Feynman-Hellmann theorem does not reproduce 
a correct answer for the individual contributions of the valence and Dirac-sea quarks to 
the nucleon scalar charge, even though the net result, i.e. the sum of them, is correctly 
reproduced. 



C. A careful treatment and a resolution of the puzzle 

To reveal the origin of the discrepancy above, we first recall a general proof of the 
Feynman-Hellmann theorem in the form convenient for our discussion below. We start 
with the expression of the energy given as 

E{a) = {^{a)\H{a)\^{a)). (36) 

Here, we assume that |^E'(a)) is normalized as {"^{a) \ ^(a)) = 1. For the proof of the 
theorem, the state |^'(a)) need not be an exact eigenstate of the Hamiltonian H{a). For 
instance, it can be an approximate eigenstate in Hatree-Fock theory, which is variationally 
optimized with respect to the Hamiltonian [9]. Under a small variation of a parameter a, 
the change of E{a) is given by 

6E{a) = {^{a)\6H{a)\^{a)) 

+ {■^{a)\H{a)\6^{a)) + ((5^(a) | i7(a) | ^(a)). (37) 

If the state |^E'(«)) is variationaly optimized with respect to the Hamiltonian, the 2nd line 
of the above equation is expected to vanish, i.e. 

{^{a)\H{a)\S'^{a)) + c.c. = 0, (38) 

where, c.c. means the complex conjugate of the 1st term. We therefore obtain 

6E{a) = {^{a)\6H{a)\'^{a)), (39) 

or equivalently 

— E{a) = {^{a)\—H{a)\^{a)), (40) 



which proves the celebrated Feynman-Hellmann theorem. 
What happens if the Hamiltonian consists of two terms as 

H(a) = Hi(a) + H^ia). 



(41) 



Here, we are imagining the decomposition of the total energy into the contribution of the 
valence quarks and that of the Dirac-sea quarks in the CQSM. Note that, in the CQSM, 
this decomposition can in fact be realized as follows : 



H, = {0\H\0) 4ao, 
H2 = {^\H\n) ai an, 



(42) 
(43) 



by taking the eigenstates of the Dirac Hamiltonian as a complete set. Now, the change 
of E{a) under the variation a ^ a + 5a is given as 



SE{a) = 5Ei{a) + 5^2 (a) 

= (^(a) I 5/7i(a) + 5H2{a)\^{a)) 

+ {^{a)\Hi{a) + i^2(a) 1 5^(a)) + c.c. 



(44) 



Assuming that |\l'(a;)) is variationally optimized with respect to the total Hamiltonian, it 
still holds that 

(^(a)|ifi(a) + H2{a)\5-^{a)) + c.c. = 0. (45) 
However, this does not necessarily mean that Hi term and H2 term separately vanish as 



{^){a)\Hi{a)\5^{a)) + c.c. = 0, 
(^(a) |if2(«) |5^(a)) + c.c. = 0. 



What is meant by fj45l) is only the identity : 

{^{a)\Hi{a)\5^{a)) + c.c. = - [ (^(a) | if2(a) | 5 ^(a)) + c.c. 
On account of this observation, we therefore propose a decomposition. 



, ^ , , , dH2{a) I , / N\ 



d_ 

da 



Eo(a) - 



{^{a)\Hi{a)\^^) + c.c. 
oa 

(^(a) I H2ia) I —Q^) + c.c. 



(46) 
(47) 



(48) 



(49) 
(50) 
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Applying this result of general consideration to our case of interest, we obtain 



= (Tj,al + CTsea, (51) 

with 

^val = ^ial^'' + Sa^ah (52) 

^sea = ^^JeP + ^^sea- (53) 

Here, the a^^^^ terms correspond to the answer obtained with naive application of the 
Feynman-Hellman theorem as explained in the subsection A, i.e. 

= ^^-'(^o), (54) 
= (55) 

On the other hand, the correction terms to this naive answer is given by 

{^{mo + Amo)\H\^{mo + Amo)y''^ - (mo) \ H \ ^ (mo))'"'^ 

oo^ai = - lim , (56) 

Amo-»0 Amo 

(^(mo + Amo) |//|^(mo + Amo))'"'' - (^(mo) I if I ^(mo))"^'' 

ocJsea = - ^lim , (57) 

Amo^O Amo 

with the simplified notation H = H{mQ). We emphasize that the identity fl48p dictates that 
6avai and 6asea are not independent but must satisfy the constraint : 

6ayai + Sa^ea = 0. (58) 

That is, the above correction terms generally contribute to both of the valence quark term 
and the Dirac-sea term, but they are expected to cancel in the sum, i.e. in the net contri- 
bution to the nucleon scalar charge, or equivalently in the nucleon sigma term. 

Since all the quantities appearing in the above discussion can be calculated explicitly 
within the CQSM, we can verify whether our theoretical consideration is correct or not. 
Shown in table I are the results of our numerical calculation for the relevant quantities. 
The 2nd row of the table show the contributions of the valence and Dirac-sea quarks to the 
quantity a'^^^^ together with the sum of them, while the 3rd row give the corresponding 
contributions of correction term 6a. One sees that the valence quark contribution to 6a is 
large and negative but the Dirac-sea contribution to 6a has just the same magnitude with 
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TABLE I: The CQSM predictions for the nucleon scalar charge. The calculation by using the 
Feynman-Hellmann theorem is compared with the direct calculation. 





valence 


Dirac-sea 


total 


aim 


11.18 


- 4.31 


6.87 


6a 


- 9.27 


9.25 


- 0.02 


a(^^^+6a 


1.91 


4.94 


6.85 


a(D) 


1.91 


4.95 


6.86 



opposite sign (aside from very small numerical error). The 3rd row represents the sum of 
a^^^^ term and 5a term, whereas the 4th row stands for the answer of the direct calculation 
of the nucleon scalar charge. One can clearly convince that, if one properly takes account of 
the correction term Sa in addition to the term a^^^^ naively expected from the Feynman- 
Hellmann theorem, the answers of the direct calculation is legitimately reproduced not only 
for the net scalar charge but also for the individual contributions of the valence and of the 
Dirac-sea quarks. 

Now, we have confirmed that, within the framework of the CQSM, the direct calcu- 
lation and the indirect calculation by utilizing the (slightly modified) Feynman-Hellmann 
theorem give exactly the same answer for the decomposition of the valence and Dirac-sea 
contributions to the nucleon sigma term. The answer clearly shows the dominance of the 
contribution of the Dirac-sea quarks over that of the valence quarks, in sharp contrast to 
the corresponding answer of the lattice QCD simulation with use of the Feynman-Hellmann 
theorem The lattice QCD version of the Feynman-Hellmann theorem is derived based 

fi n 

on the framework of partially quenched QCD (PQQCD) p|,[l8[, where the valence quarks 



that coupled to external sources for the asymptotic hadrons are distinguished from the sea 

□ n 

quarks that contribute to the quark determinant [20|,[21|. By treating the masses of the 
valence and sea quarks as independent variables, the PQQCD version of Feynman-Hellmann 
theorem is written down in the following form : 

1^ = {N\UU + dd\N),onn, (59) 
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= {N\uu + dd\N)disc, (60) 

where the short-hand notation to omit the vacuum subtraction term — V{0\uu + dd\0) is used 
for the disconnected piece. A general strategy for evaluating these terms are as follows. One 
first generates statistically independent ensembles of gauge field configurations at several 
different sea quarks masses. After that, one measures the nucleon mass for various valence 
quark masses on each of those gauge ensembles. This in principle makes it possible to 
evaluate valence and sea quark mass dependence of the nucleon mass. The physical answer 
for the nucleon sigma term is then obtained by calculating the derivatives fl59|) and fl60|) of 
M^v at the unitary point rrisea = ^vai- (In practice, the simulation in the chiral region is not 
economical, so that the results of simulations in the larger quark mass region are extrapolated 
to obtain answers corresponding to the chiral region with the help of the partially quenched 
baryon chiral perturbation theory [22].) 

It appears that there is no question about this general prescription. How can we reconcile 
the prediction of the lattice QCD with that of the CQSM, then ? Naturally, an easy 
explanation is to claim that the decomposition of the valence and Dirac-sea contributions to 
the nucleon sigma term in the CQSM does not simply correspond to that of the connected 
and disconnected contributions to the same quantity in the lattice QCD. We cannot deny 
this possibility completely, because there is no rigorous correspondence between the two 
theories and their decompositions of the nucleon sigma term. From a physical viewpoint, 
however, the discrepancy seems too large (or more than quantitative) to accept this naive 
conclusion. In our opinion, this discrepancy should be taken more seriously. If there is any 
resolution to this problem, we conjecture that it must be traced back to a difference between 
the treatment of the valence and sea quarks in the CQSM and that in the partially quenched 
QCD. In our framework of the CQSM, we do not need to distinguish the masses of the valence 
and Dirac-sea quarks. They are treated on the equal footing from the beginning to the end. 
On the other hand, there is an apparent asymmetry in the treatment of the valence and 
sea quarks in the PQQCD. (We are talking about the asymmetry in the treatment of the 
valence and sea quarks when generating ensemble of gauge field configuration. Ideally, all 
the field configurations of the nucleon constituents, i.e. the gluon field, the valence and 
sea quarks, should be determined according to a self-consistent dynamics of QCD.) It may 
be certainly true that the masses of the valence and sea quarks are taken equal at the 
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end of calculation, and that the physical answers obtained in this unitary limit is taken to 
be physical [2I]. This would also apply to the PQQCD version of the Feynman-Hellmann 
theorem. In consideration of fundamental importance of the problem, however, we think 
it very important to check the validity of the Feynman-Hellmann theorem in an explicit 
manner by carrying out a direct calculation of the connected- and disconnected-diagram 
contributions to the nucleon sigma term within the same framework of the overlap fermion. 
(The direct calculation means a calculation of the three point functions with an insertion of 
the scalar density operators.) 

Before ending this subsection, it may be useful to recall one plausible argument, which 
strongly indicates that the contribution of the valence quarks cannot be a dominant term 
of the nucleon sigma term |23|]. As is well known, the recent analyses of the pion-nucleon 
scattering amplitude favor fairly large nucleon sigma term ranging from 50 MeV to 70 MeV 



2J]-[26|. Depending on the uncertainty of the average u- and (i-quark masses, this implies 
fairly large nucleon scalar charge a of the order of 10. As we shall argue below, it is unlikely 
that such a large value of a can be explained by the contribution of three valence quarks 
alone. To convince it, let us consider a relativistic bound state oi Nc{ = 3) quarks. Assume 
that these quarks are confined in some mean field or confining potential. A typical example 
is the famous MIT bag model. The ground state wave function of this popular model is 
given as 

%Ar) ^ [ ^<^>^- ] . (61) 

where /(r) and g{r) are the radial wave functions of the upper and lower components, while 
Xs is an appropriate spin wave function. The nucleon scalar charge in this model is easily 
obtained as 

a = {N\uu + dd\N) = N, \{f{r)f - ig{r)Y] dr, (62) 

JO L -I 

with R the bag radius. Undoubtedly, the magnitude of this quantity is smaller than Nc, 
since the radial functions satisfy the normalization, 

r + i9{r)f] r'dr = 1. (63) 

It is clear that this observation does not depend on the exact form of the mean field or 
the confining potential, so that it is quite general. As a consequence, for any model of the 
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nucleoli, which contains Nc valence quark degrees of freedom alone, we must conclude that 
there exists a upper bound such that 



ayai < Nc. 



(64) 



Our prediction in the CQSM, i.e. a^ai — 1-91, as well as the direct calculation in the 
qeunched lattice QCD in [lO|, i.e. {N \ uu + dd \ N)connected — 2.323(15), satisfy the above 
bound. On the other hand, the recent estimate by the JLQCD collaboration utilizing the 



Feynman-Hellmann theorem [ij, i.e. {N\uu + d d \ N) connected — 5.27(75) — 7.92(8), lies 
outside this bound. Again, highly desirable is a direct calculation of the nucleon sigma 
term within the framework of the overlap fermion, without utilizing the Feynman-Hellmann 
theorem. 



III. SUMMARY AND CONCLUSION 



In summary, we have investigated the nucleon sigma term or the nucleon scalar charge 
within a simple effective model of QCD, i.e. the chiral quark soliton model. It was demon- 
strated that the naive application of the Feynman-Hellmann theorem does not reproduce 
the correct answer for the separate contributions of the valence and Dirac-sea quarks to the 
nucleon sigma term, which can be obtained by the direct calculation within the same model. 
It was also shown that a careful inspection of the derivation of the Feynman-Hellmann the- 
orem indicates the necessity of a correction term, which fills up the gap between the direct 
calculation and the naive application of the Feynman-Hellmann theorem. Anyhow, by using 
two completely independent methods of calculation, we have confirmed that the contribution 
of the Dirac-sea quarks dominates over that of the valence quarks in this unique observable 
of the nucleon. 

This observation however appears to contradict the corresponding answer of the recent 
lattice QCD simulation by JLQCD collaboration based on the action of overlap fermion. 
They estimated the separate contributions of the connected and disconnected diagrams to 
the nucleon sigma term by utilizing the lattice QCD version of the Feynman-Hellmann theo- 
rem, which is derived within the scheme of PQQCD, and found that the connected- diagram 
gives a dominant contribution to the nucleon sigma term and the disconnected-diagram con- 
tribution is of secondary importance. Although we do not have any convincing reasoning 
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to suspect the validity of the lattice QCD version of the Feynman-Hellmann theorem, it is 
highly desirable to check the vahdity of it by a direct calculation of the nucleon sigma term 
within the same framework of overlap fermion. This is especially so, because the separation 
of the nucleon sigma term into the contributions of valence and sea quarks seems to be a very 
delicate operation as our model analysis has shown, and also because the direct confirmation 
of the theorem is of fundamental importance to check whether the theoretical framework 
of the PQQCD, which was invented for handling loops of sea quarks in the lattice QCD, is 
working as it is expected. 
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